Coordinate Representation of the Two-Spinon wavefunction and Spinon Interaction 
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By deriving and studying the coordinate representation for the two-spinon wavefunction, we show 
that spinon excitations in the Haldane-Shastry model interact. The interaction is given by a short- 
range attraction and causes a resonant enhancement in the two-spinon wavefunction at short sep- 
arations between the spinons. We express the spin susceptibility for a finite lattice in terms of the 
resonant enhancement, given by the two-spinon wavefunction at zero separation. In the thermo- 
dynamic limit, the spinon attraction turns into the square-root divergence in the dynamical spin 
susceptibility. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.50.Ee 
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I. INTRODUCTION 

One of the most important issues in contemporary 
physics is spin fractionalization, which takes place in 
strongly-interacting one dimensional (ID) antiferromag- 
nets. At odds with what one would expect, the elemen- 
tary excitation of these systems is not a spin-1 spin wave, 
but a gapless spin- 1/2 excitation: the spinon. 

Since their discovery within the framework of the 
Heisenberg model (HM), by Fadeev and Takhtajan |Q, 
spinons are believed to be noninteracting particles since, 
according to the Bethe-Ansatz solution of the HM 
the energy of a many-spinon solution is apparently given 
by the sum of the energies of each isolated spinon. In 
this paper we challenge the non-interacting spinons idea 
through a careful analysis of the two-spinon dynamics 
in an exact solution, which shows that they do actually 
interact. 

Extensive amount of theoretical studies proved that 
spin fractionalization is a generic phenomenon in one- 
dimensional spin- 1/2 interacting antiferromagnets ||. 
Indeed, the large-scale physics of all these systems is 
always the same, given by "spinon gas" dynamics |Q. 
Therefore, without any loss of generality, one can choose 
to analyze a particular model, where the excitations are 
easier to visualize. 

The simplest model that describes the properties of 
spinons is the Haldane-Shastry model (HSM) ||[| , where 
spins- 1/2 located at the sites of a circular lattice antifer- 
romagnetically interact and the interaction is inversely 
proportional to the square of the chord between the two 
sites. In this paper we investigate the basic features of 
the HSM by employing a formalism based on analytic 
variables on the unit radius circle 0. By using real 
space coordinates, the spin-1/2 excitations become easier 
to construct and visualize than by making use of plane 
waves [|]]|]. The formalism can be easily generalized to 
the study of other models. Such a formalism allows us 



to write a "real space" representation of the two-spinon 
wavefunction. 

By analyzing the real space two-spinon wavefunction, 
we show that spinons scatter by means of a short-ranged 
attractive potential and analyze in detail the physical 
consequences of the existence of this potential. The 
short-rangeness of the interaction makes spinons free 
when they are widely separated. However, from the exact 
solution of the Schrodinger equation for two spinons we 
find that the amplitude of the wavefunction is greatly en- 
hanced when they are on top of each other, phenomenon 
which we refer to as "resonant enhancement" . While the 
density of states is uniform at low energy, resonant en- 
hancement causes the overlap between the wavefunction 
for the localized spin wave and that for the spinon pair 
to be significant, but not enough to create a two-spinon 
bound state. The corresponding matrix element is en- 
hanced so as to make the spin-1 excitation absolutely 
unstable. 

Physical consequences of the instability of the spin- 
wave appear in the functional form of the dynamical 
spin susceptibility (DSS), X?( w )- The DSS is the Fourier 
transform of the spin-spin correlation function. Its func- 
tional form can be experimentally tested by means, for 
instance, of neutron scattering experiments, the probed 
quantity being the spectral density of states, l/7rlmx g (w) 
|l0f . A system with a stable spin-1 excitation would show 
a sharp pole in Imx g (a;) at the corresponding dispersion 
relation, u> = u>(q). On the other hand, absolute insta- 
bility of the spin wave against decay into spinons will 
generate a branch cut in lmx{q,uj) at the threshold en- 
ergy for the creation of a spinon pair, which is a signal of 
the opening of a decay channel, corresponding to the lack 
of spin wave integrity. Consequently, a sharp square-root 
singularity shows up at the threshold for the creation of 
a spinon pair, on top of the broadening in the spectral 
density of states. Experiments performed onto quasi one- 
dimensional antiferromagnets provide clear evidence for 
broad spectra, while no sharp spin-1 resonance has been 
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seen JTo| ] . 

An exact calculation of the DSS cannot, in general, 
be performed, even for models exactly solvable with the 
Bethe-Ansatz. However, the HSM has the remarkable 
property that the wavefunction for a spin-1 excitation 
is fully decomposed in the basis of the two-spinon eigen- 
states . This allows us to write an exact expression for 
the DSS even for a finite lattice, thus letting us explicitly 
show the relationship between the resonant enhancement 
(Fig.(g)) and the DSS. 

The paper is organized as follows: In section II we 
shortly review the HS Hamiltonian and its exceptional 
symmetry; In Section III we introduce the ground state 
of the HSM and its representation as a function of an- 
alytic variables on the unit circle. In terms of the an- 
alytic variables the ground state takes the same func- 
tional form as the fractional quantum Hall wavefunction, 
which corresponds to a nondegenerate disordered spin 
singlet. We discuss at length several properties of the 
ground state, how to derive the corresponding energy 
and the meaning of the disorder in the ground state; In 
section IV we analyze the one-spinon solution and de- 
rive its relevant properties; In section V we focus onto 
the two spinon solution. We derive the energy eigenval- 
ues, the corresponding eigenvectors, and their norm. A 
discussion about spinon statistics is provided at the end 
of the Section; Sections VI and VII contain the key re- 
sults of our work; In Section VI we write the Schrodinger 
equation for the two-spinon wavefunction, whose solu- 
tions are hypergeometric polynomials. From the behav- 
ior of the two-spinon wavefunction, we infer the nature 
of the interaction between spinons: a short-range attrac- 
tion. The physical consequences of such an interaction 
are discussed at length in Section VII, where we derive 
an exact closed-form expression for the dynamical Spin 
Susceptibility in terms of the two-spinon wavefunctions 
and rigorously prove that the DSS is fully determined by 
spinon interaction. In the thermodynamic limit spinon 
interaction turns into the square root divergence in the 
DSS; In Section VIII we provide our main conclusions. 

II. HALDANE-SHASTRY HAMILTONIAN 

The Haldane-Shastry model |5||| is defined on a lattice 
with periodic boundary conditions. Let TV be the num- 
ber of sites. Let z a , with = 1, be a complex number 
representing a lattice site on which a spin-1/2 electron 
resides, and let S a be a Heisenberg spin operator act- 
ing on that electron. The Haldane-Shastry Hamiltonian 
takes the form: 

\ N ' ^ Za ~ z ^ 

where J is the coupling strength. The interaction is an 



analytic function of the coordinates. This is related to 
the property of a complex variable z laying on the unit 
circle: z* — z^ 1 which implies: 

1 _ ZqZ0 

\z a -zp\ 2 {z a -zp) 2 

The representation in terms of the analytic variables z a , 
which we will use throughout the paper, comes out to 
be very useful for describing the properties of spinons in 
real space. The Hamiltonian in Eq. (|l|) is clearly invariant 
under spin rotations generated by the total spin: 

N 

[H HS ,S]=0 S = ^S a . (2) 

a 

It also possesses an additional symmetry generated by a 
vector operator independent on S: 

[H„ s , A]=0 A=ij2(^ ± ^)(S a xSp) . (3) 

1 a ^ Za Z P J 

That A commutes with Tins can be seen as follows: 
[HhsA 

= £ £ V-z,\z -z,p ^ x ^ & ' ^ 

-4 ? V z _l±It 1 [(£ . § k ) Si - (S e ■ S k ) Sj 

+*E ?!! fc u L 12 . (4) 

j^k ■? Zk ' Zj Zk ' 

Although they both commute with Tins, S and A do not 
commute between themselves, being that A is a vector, 
as shown by the commutation relations: 

[S a , S b ] = i e abc S c [S a ,A b ] = i e abc A c , (5) 

From the commutators in Eq.(||) it follows that T-Lhs, 
S 2 , and (A • S) all commute with each other. The ex- 
tra symmetry of TLhs is the reason for the exceptional 
degeneracy of the energy eigenstates, as pointed out and 
discussed in [pi. The algebra generated by the two vec- 
tor symmetries of Tins is referred to as Yangian and is 
discussed in . A can be physically interpreted as the 
spin current operator for the HSM, as we show in Ap- 
pendix C. 

Starting from next section we will review the proper- 
ties of the ground state and of the one- and two-spinon 
excited states of the HSM. This will allow us to to de- 
fine the formalism we will use in order to describe the 
relevant physical properties of the model. 
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III. GROUND STATE 



B. Singlet Sum Rule 



Let N be even. We proceed by first giving the rep- 
resentation of the ground state |^gs) in terms of the z- 
coordinates and then proving that it is the actual ground 
state of TLhs- Y&gs) is defined in terms of its projec- 
tion onto the set of states with M = N/2 spins up and 
the remaining spins down. If z\, . . . , zm are the coordi- 
nates of the up spins, one defines the state \zi, . . . , Zm) 

as: \zi, z M ) = H j=1 Sf n a =i c lj°) wnere 1°) is tlic 
empty state. The projections are given by: 



®Gs(zi, zm) 



M 

life 

j<k 



A I 



Zk 



) 2 n 



(6) 



where Z\,...,Zm denote the locations of the | sites all 
others being |. We can imagine the spin system as 
a 1-dimensional string of boxes populated by hard-core 
bosons, the J. spin state corresponding to an empty box 
and the | spin state corresponding to an occupied one. 
The total number of bosons is conserved, as it is physi- 
cally the same thing as the eigenvalue of S z . Let us, now, 
review the main properties of ^gs(zi, . . . , zm)- 



A. The norm of ^gs 

^Gs(zi, . . . , zm) is a homogeneous polynomial of de- 
gree N — 1 in the variables z\,..., Zm- Its norm can be 
computed by using the following identity: 



C 



M 



M 

E 

XI,. ..,2m l<j 



= (£)«/^,.i!^TT(i-^ , ( 7 ) 
2m J zi J z M A , A Zj 



where the integrals are calculated on the circle of radius 
1. The integral in Eq.(Q) has been evaluated by Wilson 
lOl. The result is: 



We shall prove that the ground state is a spin singlet 
by showing that |^gs) is annihilated by both S z and S~. 
S z \^gs) = because I^gs) has an equal number of | 
and I spins 



[S *GS](Z2, ■ ■ -,Z M ) = 



N 



{Z2, 



,ZM\S a |*gs) 



= g*a T,jAT, z i\^os(z 1 ,...,z M ) = (10) 

(.=1 ' ^a=l ' 1 

since Y,a=i z i = N 6 to (mod N) . 



C. Coordinate Invariance 

Spin rotational invariance implies that ^as is invariant 
under interchange of | and j coordinates. More generally, 
the quantization axis can be taken to be an arbitrary di- 
rection in spin space. Denoting the sites complementary 
to zi, z M by rji, ...,»)m, so that 



M 



Y\_(z - z k )(z - Tjk) = z N - 1 , 



we have for fixed j 

M 



\\{z } - z k )(zj - r] k ) = lim 



N z 



N-l 



and thus 



z^Zj Z — Z, 



(11) 



(12) 



M 



M 



M 



IK* 



Zk) 



j<k 



J 3,k J 



M M 

(-i) M n^-^) 2 ii^ 

j<k j 



(13) 



J_ M fd^ JdzMjr {1 _z i) 
y 2m> J z x J z M l} y z/ 



therefore: 



ft, = 



2 M 



(2M)\ 

oM 



(8) 



(9) 



D. Reality 

The ground state is its own complex conjugate, and 
therefore is a real number: 



M 



y* GS ( Zl ,...,z M ) = JJ( 

j<k 



*\2 
Zj — Z k ) 



M 

n 



M 

life 

j<k 



M 



''\{z)- N =*gs(zi,...,z m ) ■ (14) 
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E. Translational Invariance 



G. Ground State Energy 



The crystal momentum of the state, q, is defined (mod 
2tt) by the equation: 



^Gs(ziz, . . . , z M z) = e lq ^ GS {zi, . . .,z M ) 



(15) 



where z=exp(27ri / N) . From Eq. (|15|) it comes out that q 
can be either or it, according to whether N is divisible 
by four or not. ^gs equals itself, up to an overall minus 
sign, when translated by one lattice constant. 



F. Disordered State 

|^Gs) is a disordered state. The way spin-spin cor- 
relations fall off with the distance defines whether a 
state of a magnetic system takes order or not. The 
relevant quantity is the spin-spin correlation function, 
X{z a ) = (^gs\S^S-\^ gs )/(^gs\^gs), which can be 
expressed in terms of two-spinon wavefunctions only, as 
we show in Section VII. 



0.2 p 

0- 



-0.2 



-0.4- 



tt/2 



FIG. 1. Spin-spin correlation function decay for N = 60. 

One-dimensional systems do not break continuous sym- 
metries, so they are not alleged to order. However, there 
is a substantial difference between half-integer spin chains 
and integer spin ones H . Both have a disordered ground 
state, but the former have excitations above the ground 
state which are gapless in the thermodynamic limit while 
the latter have a gap that survives the thermodynamic 
limit and is given by A = hv/£, where v is the spin- wave 
velocity of a nearby ordered state and £ is the correla- 
tion length. The consequence of this is that the fall-off of 
the spin correlations in the ground state of half-odd spin 
chains is not as abrupt as for integer-spin chains, where 
the correlations are suppressed within one or two lattice 
spacings ("Haldane's conjecture"). Fig.(|l|) shows that 
the behavior of the HSM is the one expected for half-odd 
spin chains. Correlations decay as (— l) x /x, according to 
Haldane's conjecture. 



|*gs) is an eigenstate oiH.HS with eigenvalue: 

".2\ / 5 



Hhs^gs) = -J 



24 



(16) 



We trade sums over spins on the lattice for derivative op- 
erators that are understood to act onto the analytic ex- 
tension of ^gs(zi, • ••■> zm), in which the zj's are allowed 
to take any value on the unit circle. After computing the 
derivatives, we constrain them again to lattice sites. We 
begin by observing that [S^Sg^csKzi, Zm) is iden- 
tically zero unless one of the arguments z\,...,zm equals 
z a . We have 



N 



*GS](Z1, 



ZM) 



M N 

X/ EI u - 

3=1 P& ' 3 



Z f3\ 



r*&s(Zl, ■ • • , Zj-l,Zp, Z j+1 , . . . , Zm) 



V Z ^ Z/3 - Zj ^ ] ( - \l[ ^Gs{zu---,Z M ) 

^^\^l\\z i -Zg\ 2 \ K dz i , \ Z4 
j=i e=o K p=£j 1 3 p 1 3 3 k 3 



N-2 M t+1 M f T , , 

V^V^j „ J ^ G s(z u ...,z M ) 



7 . y. m Ai{-—s, 



(17) 



The coefficients Ai are evaluated in appendix B. Their 
remarkable property is that they are zero for N > / > 2. 
Hence, Eq.@ can be rewritten as: 



M 



3=1 k 



J (N — 1)(N — 5) N- 3 2 d 



12 



2 3 dz 



1 3 2 \f *Gs(z 1 ,...,Z M ) 

'2 Zj dz 2 J\ 



M 



N(N- l)(JV-5) 2z 3 

9 



24 



Zj Zfc 



2z 2 



M 2 
Z J_ 



•in ( z i ~ z k ){zj - z m ) 4ri (zj - z k ) 2 

j^tk^rn J j^k J 



X^ GS (Zi, . . .,Z M ) 



( AT 1U | "I 

= -T-Eu -p\*asizi,...,z M ) . (18) 



4 



In Eq.(TLg) we have made use of the rule: 



I. Degeneracy 



[Z a - Zp){z a - Z 7 ) (zp - Z a )(zp - 2 7 ) 



(Z 7 - 2 Q )(Z 7 - Zp) 



1 . 



We also have: 



JV 

£ 

p^a 



S^S 



a u P 



\z a - Zp | 2 



N(N 2 - 1) 
48 



&GS](Z1, ...,Zm) 



M 

£ 



2fc | 



(19) 



The HSM ground state is not degenerate, but is nearly 
so. We already pointed out that half-odd spin magnets 
have a gapless spectrum. In the next sections we will see 
that elementary excitations above the ground state are 
spinons and that their spectrum is relativistic. In partic- 
ular, at the endpoints of the Brillouin zone, the energy 
becomes the same as the ground state energy, modulo 
corrections that are sub-leading in the thermodynamic 
limit. This means that, in principle, one can have many 
states with the same energy as the ground state which 
are distinguished from one another by their number of 
spinons. 

An example is provided by the singlet state of two 
spinons with total momentum ir. It is given by: 



X &GS(Z1, ■ ■ -,Z M ) ■ 

This completes the proof, since 



N 



2 X N' 1 ^ \ z a -zpV 

a^P 1 P 1 



N 

£ 

a^P 



QZ QZ 
J a Jp 



Z &\ 



(20) 



M 



(21) 



The wavefunction ^gs{z\i ■ ■ ■ > z m) was first introduced 
by Haldane and Shastry [||](| in analogy to the exact 
Sutherland solution of the continuum limit of the prob- 
lem jl3| . The proof that this wavefunction is the actual 
ground state of TLhs is a consequence of the factorization 
of the HS Hamiltonian, as we are going to discuss next. 



H. Factorization of Whs 

In Appendix D we prove that Tins can be written as: 

JV 



iV(jV 2 + 5) , N + lp 



48 
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The operators D a are given by: 



1 x -\ z a 

2 ^ z~ 

Pjta 



Z & 



zp 



i(S a x Sp) + Sp 



(22) 



(23) 



and they annihilate I^gs) ( see Appendix D). Eq.(|2^) 
implies that |\&gs) is the ground state of the HSM be- 
cause TIhs can be written as a constant plus nonnegative 
definite operators, and the only state satisfying the re- 
quirement of minimum energy is |Wgs)- 



$>s(z-l,...,z m ) = Ylizj-Zk) 2 



M 



Its energy is given by: 



n HS \^s) = -J 



24 



(24) 



(25) 



and is the energy of the ground state plus corrections 
that go to zero in the thermodynamic limit. 



J. Spin current 

We now show that I^gs) is an eigenstate of A 2 belong- 
ing to the 0-eigenvalue. The action of A z on the ground 
state gives: 



A z |*gs) = . 
Eq.(p6|) can be proved as follows: 

[A.*V as ](zi,...,z M ) 



(26) 



M N 

^££ 



zp 



zp 



JV— 1 

£ 

i=o 



2 * — ' * — ' ' v Zj — zr 

3=iPjtj V 3 P 



(zj - zp) 1 d f^Gs(zi,---,z M ) 
V, dz 



N N 
--(N-2) + -(N-2) 



x Wgs(zi, ■ ■ • , z M ) = , 



(27) 



where we have made use of the results of Appendix B and 
of the technique described in detail in subsection III.F. 
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IV. ONE-SPINON WAVEFUNCTION. 



Let the number of sites N be odd and let 



At odds with the naive idea that the elementary ex- 
citations for interacting magnets are integer spin states 
(spin flips), Faddev and Takhtajan |Q first conjectured 
that one-dimensional half odd spin chains exhibit exci- 
tations given by spinons carrying half-odd spin. For a 
chain with an even number of sites, the ground state is a 
disordered spin singlet but, if the number of sites is odd, 
the minimum possible value for the total spin is 1/2. In 
the thermodynamic limit, it makes no difference whether 
one begins with an odd or an even number of sites. The 
short-rangeness of the correlations in the ground state 
makes it insensitive to the boundary conditions, so, in 
the thermodynamic limit, there is no way to distinguish 
between chains with odd or even number of sites. States 
with half-odd spin are, then, alleged to appear as eigen- 
states of Tins with an odd number of spinons. In this 
section we shall present the one-spinon wavefunction and 
discuss its properties. Following Haldane || we consider 
a wavefunction of the general form 



#(zi, ...,z M ) = ®(zi, ■ 



M 

•>ZAf) Y[(Zj 
j<k 



M 



Zk 



JJzj (28) 



where Zi,...,Zm denote the position of the up spins. 
Here $ is a homogeneous symmetric polynomial of de- 
gree less than N — 2M + 2 in each variable. This latter 
condition causes 'J to be a polynomial of degree less than 
N + 1 in each of its variables Zj , and thus allows the Tay- 
lor expansion technique used for the ground state to be 
applied. Doing so, we find that 



* Q (z 1; ...,z M ) 



M 



IK* 



M 

II <* 



M 



(31) 



where M = (N — l)/2. This is a J. spin on site a sur- 
rounded by an otherwise featureless singlet sea. It is 
worth stressing that Eq.(|3l|) makes perfect sense for any 
z a on the unit circle. Nevertheless, as z a coincides with 
a lattice site, it represents a spin J, localized at the cor- 
responding site. The spin density of the corresponding 
state, plotted as a function of the spinon position, will 
be uniformly zero, as appropriate for the disordered spin 
singlet, except for an abrupt dip centered at z = z a (see 
Fig.(^)). Such a dip is what we refer to as "real space 
representation" of a spinon at z a . Hence, a spinon can 
be visualized as a local defect in an otherwise feature- 
less singlet sea. This defect behaves like a real quantum 
mechanical particle, as we will show in the following. 



0.5 



<s!> 



P 

-0.5 



a 



1 



^ ( 4M 2 -1)-|m 2 WT, . 



(29) 



provided that $ satisfies the eigenvalue equation for A 



He 



,<9 2 $ 



M 



4z 2 a$ 



dz 

3 



j j^k Zj _ Zk dzj 



2 - 







a - (3 

FIG. 2. Spin and charge profiles of the localized spinon 
|* a ) defined by Eq.([jl|). The dotted lines are a guide to the 
eye. 



N-3 



M 



9$ 

'dzj 



(30) 



By definition, ^ a is an eigenstate of S z with eigenvalue 
-1/2. In order to prove that it is a spin- 1/2 state, we need 
to show that S~ annihilates it. Indeed, per Eq.(]lO|) we 
have: 



A. One-Spinon Spin Doublet 

We look for one-spinon and two-spinon wavefunctions 
in the functional form given by Eq. (EM) . Here we analyze 
the one-spinon wavefunction. 



N 



(32) 



which proves that ty a is the spin- 1/2 component of a spin 
doublet. 



G 



B. One-Spinon Energy 

Eq.(31) corresponds to a particular choice of <I> in 
Eq.@, given by: 

M 

$(zi, . . . ,Z M ) = $a(zi, ■ ■ ■ ,Z M ) = Y\_( Z » ~ Z j) ■ ( 33 ) 

2=1 



Eq. (pop , once written for the state $ Q , takes the form: 



d 2 A- 3 



M-z n 



d_ 

' dz a 



The eigenstate of Hhs is given by: 



(34) 



1 - 

V m (zi,...,z M ) = 22(zl) m * a (zi,...,z M ) (35) 



and the energy eigenvalue is 



J.2tt. 2 ,iV-l Ml 
+ 2 ( A } ( ~ m) ^ r 



(36) 



with < m < (A - l)/2 and A = m((JV - l)/2 - to). 



or not. The loss of half of the states available for a reg- 
ular fermion is a peculiar property of spinon spectrum. 
No negative energy states appear, i.e., there is nothing 
like an "antispinon" . One can picture a spinon as either 
an electron or a hole whose charge has been pulled out 
by the interaction. According to such a picture, a spinon 
can arise either from an electron with the same spin or 
from a hole with the opposite spin, which explains the 
halving of the Brillouin zone. 




N=9 

I h«— • — ♦ — •— H 1 

N=ll 

I* * «l 1 1* • »l 

-71 -TZ/2 7C/2 71 

FIG. 3. Top: Spinon dispersion given by Eq.([Io|). Bottom: 
Allowed values of q for adjacent odd N. 

The spinon dispersion at low energies is linear in q with 
a velocity 



C. Crystal Momentum 

The state \*& m ) is a propagating J, spinon with crystal 
momentum 

Q=l N -^( m +\) (mod 2^) , (37) 
per the definition 

& m (ziZ,...,ZMz) = exp(ig) * m (zi, . . . , z M ) , (38) 
where z = exp(2m / N) . Rewriting the eigenvalue as 

mm) = {-J(^)(N+^- ±)+E,\\9 m ) , (39) 



v. 



spinon — ^ u 



(41) 



The half-band of single elementary excitations for odd A 
are the only 5=1/2 states without extra degeneracies. 
The ground state of the odd-N spin chain is 4-fold de- 
generate and is given by |^ m ) for to = and (N — l)/2 
and their f counterparts. This corresponds physically to 
a "left-over" spinon with momentum ±7r. 



D. Spin-Current 

We now study the action of A z on the state for one 
propagating spinon. Working as for the ground state one 
gets: 



we obtain the dispersion relation 



m = k 



(mod 7r) 



(40) 



plotted in Fig.(g). Note that the momenta available to 
the spinon span only the inner or outer half of the Bril- 
louin zone, depending on whether N — 1 is divisible by 4 



A- 1 



TO. >\^r 



(42) 



and the eigenvalue of A z comes out to be proportional to 
the spinon velocity 



dE(q) 2ir J ( N - 1 
— = - — \—- m 



(43) 
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Eq.(^||) is proven by first letting A z act on the state ^ c 
denned in Eq.(|l|): 



[A z * Q ](zi,...,z M ) 



M N 

^EE 



z 



AT-l 

E 



- zg) 1 



2 ^ " p V Zj - zg I ^ l\ 

3=1 0ft V J P / 1=0 



8 f ^ajzi, . . . ,Z M ) 

dz) \ Zj 



, M M 

J-M(JV-2) + 2££ 

1 3=1 ift 



M 



*E 

3=1 



*a(«l, . . . , 2m) 



After taking M = (N - l)/2 one gets: 



AT - f 9 . , 

1 Z »-q^~ f^a(Zl, . . . , 2m) 



(44) 



(45) 



which, on the basis of the states |^ m ), gives the result 
quoted in Eq.(E3). 



V. TWO-SPINON WAVEFUNCTION 

Let us now focus on the two-spinons state. Spinons 
maintain their integrity when many of them are present. 
This does not mean that spinons are noninteracting. 
They can be separated at large distances, therefore be- 
ing asymptotic states of the system. However, they also 
scatter strongly off each other by means of a short-ranged 
attractive potential. The interaction between spinons is 
not enough to create two-spinon bound states, but it gen- 
erates a peculiar "piling-up" of the relative wavefunction 
when the two spinons are on top of each other (Fig. ([!])). 
This is the reason for the huge decay amplitude for a spin 
wave into a pair of spinons, that is the spin fractionaliza- 
tion. 

In this section we derive the two-spinon eigenstates, 
their norm and the corresponding value of the spin- 
current. Moreover, we show that the appropriate statis- 
tics they obey is neither fermionic nor bosonic. They are 
semions, i.e., particles with 1/2 fractional statistics. 



A. Two-spinon Energy 

Two |-spinons can be pictured as two |-spins within 
an otherwise featureless disordered sea. The state with 
two j-spinons centered at z a and zg, respectively, is given 
by (N is even and M = N/2 - 1): 



^ a g{zi, . . . , Zm) 



M 



A I 



M 



E. The Norm 

The squared norm of f m is defined as: 

(*m|* TO )= E \^m{zi,...,Z M )\ 2 • (46) 

Zi,...,Zm 

By means of a simple algebraic procedure, we generated a 
recursion relation between (^ m \^m) an d (^m-il^m-i)- 
Such a procedure can be straightforwardly extended to 
the norm of the multiple-spinon states. We discuss it at 
length in Appendix E. The induction relation is: 



{V m \* m ) (m-|)(M-m + l) 



(* m _i|* m _i) m (M-m+\) 
This recursively gives: 



(*m|* m ) = -FT! 



T[M + l]T[m + h]T[M 



r[i]r[M + \]T[m + l]T[M - m + 1] 



(47) 



Cm 



(48) 



where Cm is the over-all constant we have introduced in 
Eq.@. 



= ]J(z a - z 3 ) {zg - Zj ) Y[ {Zj -z k ) 2 Y[z 3 . (49) 

3 j<k j 

As for the one-spinon case, z a and zg are not necessarily 
lattice sites. If they are, Eq.(p9|) represents a pair of 
spinons at z a , zg. To derive the eigenvalue equation, we 
start from a wavefunction in the form of Eq.(28), where 
we take the function $ to be equal to: 



M 



*«/3 = \\{z a - Zj)(zg - Zj) 



Eq.(pO|) can be rewritten for $0/3, yielding: 



(50) 



. M 



3 = 1 



2 d 2 <P a! _ 
' dz 2 



M 

^E 

3#fe 



4z 2 d<S> a g\ N-3 



Zk dzj 



M 

E 



2 ^ "• dzj 

3 = 1 



2 d 2 2 d 2 

Z~ „ „ Za 



Z a - Zg dz a Zg - Z a dzg a dz 2 a dz 2 



+ ( 



N — 3 . 



l dz a 



Z 



dzg 



2M — M{N — 2) 



a0 
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= A$ Q/3 . 

Let us now define the states ^ mn as follows: 

* m „(zi, . . .,2m) 



(51) 



In terms of spinon momenta, the expression 
ergy is: 



of the en- 



JV 

E 

a,0 



(4T (*g) T 

N N 



(52) 



A set of linearly independent states may be constructed 
by taking only the ^J mn with M > m > n > 0, which 
shows the overcompleteness of the set of states ^f mn . On 
such a set of states Eq.(|5l|) becomes: 



1/y 2j?_ 4 *j g 



M d 

3 a i 3 



3 J^k 



3=1 



-^(^-^F + ^) + ™(y-i-™) + "(y-i-") 



* m „-^(m-n+2£) # m+ i,n-* , (53) 

where £m = n if m + n < M, = M — m if m + n > M 
and, in deriving Eq. (|53|) we used the identity 



x + y 
x-y 



( rn n ri m\ 

{x y - x y ) 



2 x m - l y n+t - (x m y n + x n y m ) 



(54) 



e=o 



We look for solutions to Eq.(J53j) which are linear combi- 
nations of the states & m +t tn —i: 



(55) 



The coefficients ag are found to be: basis of the states 



-(m-n + 2£) 



2£(£ + m-n + l/2) 



J2*k"\; («o = l) (56) 



fc=i 



and the corresponding two-spinon energies are given by: 



+ [E(q m ) + E{q n ) - 



w J \q 



m Hn 



N 



(t?m < Qn) ■ (58) 



E mn is the sum of the ground-state contribution, Eqs = 
-J(tt 2 /24)(7V + 5/iV), and E(q m ,q n ), which is the two- 
spinon energy above the ground state. E(q m ,q n ) is the 
sum of the energies of two isolated spinons plus a negative 
interaction contribution that becomes negligibly small in 
the thermodynamic limit. 

Such a simple solution for the two-spinon problem is 
possible because the matrix to which Eq. ( |53"|) corresponds 
is lower triangular, i.e. takes the form 



Matrix 



So 

wio Si 

v 2 o v 2 i E 2 

«30 "31 ^32 E 3 



(59) 



where 

Erj = A 

and; 
A 



j - "m+3,n-j and v pq = -(m - n + 2p + 2q) 



N . N . m — n 
m{— - 1 - m) + n{— - 1 - n) — 



(60) 



The eigenvalues of such a matrix are its diagonal ele- 
ments, and the corresponding eigenvectors are generated 
by recursion. 

The transformation in Eq.(|55"|) can be inverted and it 
takes the form: 



(61) 



1=0 



The coefficients 6™ n can be expressed in a closed-form 
formula in terms of the coefficients of the two-spinon 
wavefunctions. We provide their expression in Section 
VI. 



J 2tt 2 
2 l iV J 



tt 2 1 Q 24 

E mn = -J—(N- — + 

24 V N N 2 ' 



r N . ,N . m — n 
m(— — 1 — m) + n(— — 1 — n) — 



(57) 



B. The Norm 



The squared norm of the state $ mn is defined as: 



\^mn{zi,...,z M )f 



(62) 
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As for the one-spinon wavefunction, we calculate the 
norm of the two-spinon states by means of mathematical 
induction. The details of our calculation are discussed in 
Appendix E. The basic induction relations are given by: 



-l*r. 



(n - §)(M -n+ §)(m - n + 1) 



l)(m 



)(m 



(63) 



m — l,n I ^m— 1,ti 



(M — to + 1)(to — n + \){m — n — h)m 
(M - m + |)(m - n) 2 (m + |) 



(64) 



From Eq.(63|6J) one gets the formula for the squared 
norm: 



mn I ^mn, 



= c 



M+\ r[m - 7i + |]r[m - n + §] 
7r r 2 [rn — n + 1] 



I, N 



(67) 



The phase to "exchange" the spinons by moving a all the 
way around the loop is thus 



Acj) 



1 ,„ N . dz a 7r , , n , 

2 (1 ~y ) 7^ = ± 2 l ( mod2 ^) ■ ( 68 ) 



This number is or n for bosons or fermions. The number 
of states available to t J. spinons, determined by counting 
the number of distinct symmetric polynomials of the form 

®z Al ,...,z Ai (zi, Z(jV_f)/ 2 ) 

(N-e)/2 

Yl ( Zj - z Al ) x ...x (z 3 - z Ae ) , (69) 



MJ 



N/2 + £/2. 
\ o > 



(70) 



r[m + 1]T[M - m + i] T[n + |]r[Af - n + 1] 
T[m + f ]r[M - to + 1] f> + l]r[M - n + |] 



(65) 



Eq. ( |65|) basically agrees with the result quoted in [jl4| , al- 
though we derived it by making direct use of the operator 
Ti-HS (See Appendix E). 



This is just halfway between the numbers 



A/? ermi = [ N { 2 ) 



(71) 



likewise calculated assuming that the number of states 
available for one particle is N/2. 



C. Spin Current 

The are eigenstates of A z . Indeed, a manipulation 
similar to the one-spinon case yields: 

A z |* m „) = -m-7iW mn ) , (66) 

with the eigenvalue given by the sum of the two spinon 
velocities. We will just skip the proof of Eq.(|66|) which 
works exactly like the proof of Eq.(f42]). 

D. Spinon Statistics 

Spinons are semions, i.e. particles obeying 1/2 frac- 
tional statistics. Since the 2-spinon wavefunction \E" a p 
has the property 

^* a p(z 1 ,...,z N/2 _ 1 ) = {z a Zp) 1 ~ N/2 ^ a p(zi,...,z N/2 _ 1 ) , 

the Berry phase vector potential for adiabatic motion of 
spinon a in the presence of /3 is 



VI. SCATTERING RESONANCE 

In this Section we provide one of the key results of 
our study: the analysis of the interaction between two 
spinons. First, we properly define the real space repre- 
sentation for the two-spinon relative wavefunction. Then, 
we study the behavior of the corresponding amplitude as 
a function of the spinon separation. Here we construct 
the real space wavefunction for a spinon pair and we show 
that our results provide a clear evidence for spinons being 
interacting particles. 

The real space representation for the two-spinon wave- 
functions corresponding to the energy eigenstate |$ mn ), 
z™Zpp mn (z a / Zp), is defined by the decomposition of the 
state of two localized spinons at z a and zp, l^ap), in the 
basis of the fully polarized two-spinon eigenstates: 

M m 

= £ £(-ir+"«?w(^)<IW • (72) 

m— n—Q ^ 

l^rnn) is an eigenstate oiTLus with eigenvalue E mn . This 
implies 
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{$mn\n HS \^ 



af3) — E mn {<fr mn \^> a ()) 



(73) 



From Eq.(|oT|) we see that (^ m n\Ti-Hs\^ap) can be written 
as a differential operator acting on {<& m n\^ap)- *q/3 is 
perfectly defined for any z a , zp on the unit circle, so, 
the differential operator acts on the analytic extension of 

<$mn|*Q/3) as 

inn 



1 z Q + ^ d d 
2zt,-z,3 az a dz/3 



(74) 



In the differential operator in Eq. ([74]) we recognize the 
sum of the energies of the two free spinons and a 
velocity-dependent interaction, which diverges at small 
spinon separations. Eq.(f74|) allows for determination of 
the exact expression of p m n{z a / zp). Indeed, by using 
Eqs.( [72| , [73|j7^ ), we find the following equation for p mn {z) 
(z = z/3/z a ): 



z(l 



d 2 p n 



(— m + n + -)z 



dz 2 



- - to + n 



dz 



-Pn 







(75) 



This equation is a special case of the hypergeometric 
equation |15| where the parameters c, b,a are given by 



1 

C = TO 

2 



b = — a = —m 
2 



The solution is a hypergeometric series whose regular so- 
lution stops at a power of z given by z m ~ n , thus becoming 
the "hypergeometric polynomial" 



Pmn (z) 



r[m — 11 + 1] 

r[|]r[m - n + |] 



m—n 



T[k- 



Mm- 



^ o r[fc + i]r[TO-n- 



fc + i] 



(76) 



The value of the spinon wavefunction at zero separation 
between spinons can be computed by means of general 
identities among hypergeometric series [ fL5| . It is given 
by: 

Pm»(l) = r[l/2]r[m - n + l]/r[m - n + 1/2] . (77) 



According to Eq.(76), |p m n(^)| 2 is the density of prob- 
ability for two spinons as a function of the distance be- 
tween them. 
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vr/2 



tt/2 

e e 

FIG. 4. Left panel: square modulus of the two-spinon 
wavefunction as a function of the separation between the two 
spinons for N = 1000. Right panel: the same plot on a log-log 
scale. The dashed line is a guide to the eye. It is a plot of 
1/x with an appropriate offset (z — exp(i#)). 

The interpretation of our results is straightforward. 
Spinons do actually behave like real particles. Indeed, 
we have been able to determine a differential equation for 
Pmni which for the two-spinon wavefunction is the same 
as the Schrodinger equation for a pair of ordinary parti- 
cles. The interaction between spinons is clearly shown in 
Fig.(||), where we plot \p m n{z)\ 2 ■ At large separations, 
the probability density oscillates and averages to 1, in- 
dependently on the distance between the spinons. This 
is a typical feature of nonintcracting particles. Indeed, 
at large separations spinons are in fact noninteracting. 
However, at small separations, |p m „| 2 shows a resonant 
enhancement, which corresponds to a huge increase of 
the probability of configurations with the two spinons on 
top one of each other. The enhancement is highest at 
relative spinon momentum of n. Such features are clear 
evidence of the interaction among spinons, which can be 
characterized as follows: 

1. It is attractive: it favors configurations with 
spinons on top one of each other; 

2. It is short ranged: spinons are free at large dis- 
tances. 

As N gets larger, the resonant enhancement peaks up 
fill . Hence, the resonant enhancement safely survives 
the thermodynamic limit, even though, in this limit, the 
energy for the two-spinon solution is the sum of the ener- 
gies of the two isolated spinons. However, the attraction 
is not strong enough to create a two spinon bound state, 
even in the thermodynamic limit. This corresponds to 
the absence of a low-energy stable spin-one excitation, 
which is a typical feature of ID spin-1/2 antiferromag- 
nets. 

The attractive force may also be inferred from the en- 
ergy eigenvalue if we rewrite it as 



E n 



E, 



Qn 



v„. 



-In 



|4>, 
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-J(^)(N+^) + E(q m ,q n ))<l> 



where 



Va 



-4^ 



(78) 



(79) 



Note that this potential vanishes as N — > oo, as expected 
for particles that interact only when they are close to- 
gether. However, as we already pointed out, the van- 
ishing of the interaction potential in the thermodynamic 
limit does not mean that no residual effects survive such a 
limit. The resonant enhancement when the two spinons 
are on the same site does survive the thermodynamic 
limit and it is the main reason of the instability of the 
spin-1 spin- wave, as we discuss at length in the next sec- 
tion. 

Before concluding this Section, we provide the expres- 



sion of the coefficients 6™ n in Eq. (|6l|) . From Eq.((j 
is straightforward to prove that: 



it 



bT n = 



T[m-n + 2£+l\ 



T[t+\]T[m- 



^ 



T[i + l] r[m- n + £+ 1] 



VII. SPIN SUSCEPTIBILITY 



(80) 



s: 



(z* a ) k (S x a - iSl) (g = 2nk/N) 



(81) 



A peculiar property of the HSM is that a spin flip at 
z a is the same as a spinon pair at the same site [^TJ. 
Therefore, we can fully decompose S~^cs in the basis 
of the two-spinon eigenstates: 

S-^ G s = ^(O^aa 



M m 

= NJ2 ^- l ) m+n Prnn{l)6{m + n- k) $„ 

m— n— 

The susceptibility is given by 

\{X\S-\* GS )\ 2 



(82) 



X 



(X\X) {^ G S\^GS) 



2(Ex — Eqs) 



{u + ir)) 2 -{E x -E gs ) 2 



(83) 



( \X) is an exact eigenstate of Hhs with energy Ex) 
Then, from Eqs.( |82] , p3| ) we have that x<?( w ) takes a 
nonzero contribution only if \X) = |$ mn ). Then, Eq. 
becomes: 



M m 

2 V V 

7U—0 n—0 v 1 



■PmJX) 



In this Section, we work out the dynamical spin sus- 
ceptibility for the HSM. We show that the DSS depends 
only on the p mn 's calculated at z = 1, which allows us to 
obtain for any finite N a simple closed-form expression 
for the DSS and to relate it to the spinon interaction. 
By carefully taking the thermodynamic limit of our re- 
sult, we obtain Haldanc-Zirnbauer formula for the DSS 
in the thermodynamic limit [ ftl| . Haldane-Zirnbauer for- 
mula shows that there is no low-energy spin-1 pole in the 
DSS, but the function takes a sharp square-root singu- 
larity at the two-spinon threshold on top of a branch cut, 
corresponding to the lack of integrity of the spin-one ex- 
citation. Our analysis definitely proves that the square 
root sharp edge on top of the broad spectrum is noth- 
ing but the interaction between spinons. The resonant 
enhancement is the square root singularity in the spin- 
susceptibility. This result is of the outermost importance, 
since it represents a way to experimentally test interac- 
tion among spinons in ID. We will come back to such a 
point in the concluding remarks. 

Let us begin with the calculation of the spin suscepti- 
bility for a finite lattice. The DSS is the dynamical prop- 
agator for a spin-1 spin flip. A spin flip with momentum 
q is created by acting on |^gs) with S~, defined as: 



x 5(m + n — k) 



2(E mn - E GS ) 



{uj + try) 2 - (E mn - E GS ) 2 



(84) 



Eq.(p4[) is another relevant result of our work. It shows 
that only the Pmn{z)'s at z = 1 determine x<?( w )- There- 
fore, the spin susceptibility is completely determined by 
spinon interaction. 

Let us analyze the thermodynamic limit of eq.(|S4|). In 
the thermodynamic limit the gamma functions can be 
approximated by using Stirling's formula: 

T[z] « V^(z- l)(*-7) e - (2 - 1) , (85) 
From Eqs.( |77| , |S5| ) we get, in the thermodynamic limit: 



N 



(*gs\*gs, 



p 2 

1 inn 



(i) 



irN(m -n+\) 



^Jn(M - n){m + \)(M - m - \) 



(86) 



Since the joint two-spinon density of states is flat, the 
sums over m and n become integrals over the (halved) 
one-spinon Brillouin zone: 
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E 



-Ml * 



A/ 



(87) 



From Eq.(||g||), we see that Eq.Q turns, in the 
thermodynamic limit, into the Haldane-Zirnbauer for- 
mula for the DSS lO 



XgH = - 



, f 91 , \Ql ~ ?2|%1 +92-?) 



y/E( qi )E(q 2 ) 



2E(q u q 2 ) 



(u + ii]) 2 - E 2 (qi,q 2 ) 



(88) 



where E(q) and E(qi, q 2 ) are the one-spinon and the two- 
spinon energies, respectively. Integration pver q\ 1 q 2 in 
Eq.([S8l) provides: 

( \ J 



x QMg) - 9 [a; - o;-i(g)] 9 [a; - ^+i(g)] , gg . 
a/w - cj_i(g)^w - w+i(g) 

and w+i(g) are the threshold energies for a spinon 
pair with momentum q, according to whether < q < tt 
or 7r < q < 2tt, respectively. They are given by 
u-i(q) = (J/2)q(n - q), tu +1 (q) = (J/2)(2tt - g)(g - tt). 
^2(9) = (J/2)q(2-K — q) is the upper threshold for the 
spin-1 excitation. From Eq.(|89|) we see that the resonant 
enhancement, given by p^ n (l), has turned into a square- 
root singularity in Xo( w ) vs - w & t fixed g, with the branch 
cut originating either at w_i(g) or at a>+i(g), depending 
on the value of q. Because the two-spinon joint density of 
states is uniform, the main conclusion we trace out from 
our calculation is that the branch cut in % g (o;), i.e., the 
broadness of the spectral density of states, is the spinon 
interaction. 

A measurement of Xqi 1 ^) m one-dimensional spin-1/2 
antifcrromagnets can be performed my means of neutron 
scattering experiments [0 . The result of the measure- 
ments [|10| does, in fact, show a sharp threshold followed 
by a broad spectrum, in good agreement with predic- 
tions of Eq. (|89|) . In light of our present discussion, we 
conclude that what is actually seen in such an experi- 
ment is a direct consequence of the spinon interaction in 
1-dimensional antifcrromagnets. Hence, the experiments 
provide evidence that spinons do interact and that the 
spinon interaction is what determines the peculiar low- 
energy physics of spin-1/2 antiferromagnetic chains. 

From Eq.(|8^) we also derive the formula for the spin- 
spin correlation function, x( z a), in terms of the two- 
spinon wavefunction at z = 1: 



(* GS \S+S;J* 

(^GS\^GS) 



Eq.(|90|) is the formula we have plotted in Fig.(|l|) 



(90) 



VIII. CONCLUSIONS 

In this paper we developed a simple approach to the 
study of spinon excitations of the Haldane Shastry model, 
based on the formalism of the analytic variables. Within 
out approach we picture spinons as local defects in the 
disordered sea. Our formalism allows for a consistent 
real-space representation of the wavefunction for two 
spinons. We construct the Schrodinger equation, whose 
solution is the two-spinon wavefunction, which shows 
that spinons behave as real quantum-mechanical parti- 
cles. By means of a careful study of the real-space two- 
spinon wavefunction, we reveal the main result: the ex- 
istence of spinon interaction and its survival in the ther- 
modynamic limit. Spinon interaction is a short-range 
attraction, which generates a resonant enhancement of 
the probability for two spinons to be at the same site. 
Such an interaction determines the low-energy physics of 
ID interacting antifcrromagnets. Since the low-energy 
joint density of states is uniform, the broadness in the 
spectral density, is exclusively caused by the resonant en- 
hancement, as we show from the finite- N expression for 
the spin susceptibility (Section VII). In the thermody- 
namic limit the resonant enhancement develops a square 
root singularity followed by a branch cut, which is the 
broadness in the spectral density of states. The branch 
cut reflects the absolute instability of the spin wave to- 
wards decay into a spinon pair. Then, we show that, 
even though in the thermodynamic limit the interaction 
is irrelevant, its main effect, the resonant enhancement, 
peaks up. 

In conclusion, we analyzed spinon interaction in an 
exact solution of the Haldane-Shastry model and its con- 
sequences for the low-energy physics of 1-D spin-1/2 an- 
tifcrromagnets. 
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APPENDIX A: FOURIER SUMS 

In this Appendix we will prove some of the formulas 
we used throughout the paper. Since the lattice sites z a 
are roots of unity we have 



N 



l[(z - z a ) = z N - 1 

a 

Then for < m < N we have 



(Al) 




C 



FIG. 5. Contours used in Eqs(|AJ) and flA3D 



JV-1 



. Z n - 1 



APPENDIX B: CALCULATIONS OF THE 
COEFFICIENTS Al. 

In this Appendix we work out the coefficients Ai that 
appear in the eigenvalue equations for the eigenfunctions 
of the HSM are defined as: 



JV-1 

E 



^ {Za - 1) 

a— 1 v ' 



2-1 



(Bl) 



They can be computed by using the equations from Ap- 
pendix A. In particular we have: 



N-l 



^o = -£ 



\ - z c 



a— 1 1 1 



(N-l)(N-5) 
12 



z a - 1 

a— 1 



a=l 



N 



z m - x dz 



N 



2wi J c (z - l)(z N - 1) 2ni J c , (z - l)(z N - 1) 



z" 1 - 1 dz 



N-l 



An 



N 
'2~Tri 



,m — 1. 
l + ( 1 )x- 



) x 2 



,N S ,Ns ,AT. , 
( x ) +( 2 )^ + ( 3 )a; 2 + 



3- 



iV-l 



A l = ~ E 4(*a - !)'~ 2 = ° > 2 ) 



(B2) 



JV+1 



(A2) 



and 



JV-l 

E 



2a -1 



JV-1 



N 



dz 



2ni Jc (z~l) 2 {z N - 1) 



N 
2nd 



dz 



c , (z-i)*(z»-i) 



1 

~2~Ki 



l+( . )x+( 



) x 2 



,N 



,N 



N 2 



( j ) + ( 2 ) x + ( 3 ) x2 + 



1 m(N — 1) m(m — 1) 



dx 
x 3 
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(A3) 



APPENDIX C: THE SPIN CURRENT 

In this Appendix we provide the physical interpreta- 
tion of the operator A, as the spin-current operator. In 
order to do so, we first construct the continuous interpo- 
lation of the lattice spin field, given by the spin density 
operator p(z). Then, we define a current density on the 
unit circle, j(z). We prove that p and j obey an equa- 
tion which, once restricted to the lattice, takes the form 
of the continuity equation for the spin density. The oper- 
ator A comes out to be the global operator whose density 
is given by j(z). 

The first step of such a construction is defining the 
field interpolating the spin operators into the interstices 
by means of the formula 



a(z) = 



y N/2 _ -N/2 



2N 



N 

E 





z + Zp 

Z- Z/3 



(CI) 



Then we can associate to a(z) a "cr model-like" Hamil- 
tonian given by: 
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1 

2iri 



da. 

dz 



da, dz 

dz z 



N 

E 



S a ■ Sp 



N " |i„ 



2/3 1 



O C2 



(C2) 



where the integral is performed over the unit circle. 
Eq.(C2) gives the HS Hamiltonian plus an irrelevant con- 
stant and an operator which commutes with it. 

We also have spin density and spin current density op- 
erators 

p(z) = —ia{z) x a(z) 



3(z) 



. a x \z^-] 
2i 1 dz 1 



^ X a 



(C3) 



That p(z) is an appropriate definition of the spin density 
may be seen by taking the limit z —> z a , being z a a site 
on the lattice. One gets 



lim p(z) = S a 



(C4) 



That j(z) is a proper spin current can be inferred from 
the continuity equation 



lim 



dj 



N 



■5 + E 



a^P 



S a ■ Sp 

Za - Zp\ 



(C5) 



The zero-momentum component of this conserved cur- 
rent density is 



1 f -. dz t 

7T- f J — = A 
2m I z 



(C6) 



The operator A is then a scaled spin current. Its action 
on the state with a fixed number of propagating spinons, 
Eqs.( p7tp2|j6^ ), is definitely consistent with such an in- 
terpretation. 



APPENDIX D: FACTORIZ ABILITY OF Hhs 



annihilates |^gs)- Indeed, by using the technique we 
developed in Section III, we find 

[Cla^GsKzi, ■ ■ - ,Zm) 



N-2 , „ N 



E 



1 ^ ZgZp{Zp - Zg) 1 \ d l f ^GS(Zg, . . . ,Z N/2 ) 



o- il p% Xa - Z f> > dz « 



N 

\ - z c . 



P±a 



Zf3 



-^ + S a ){- + S z p ) + 



N - 1 



x ^gs(zi, ■ ■ .,z M ) = 



(D2) 



for all a. However since |\&gs) is a spin singlet the irre- 
ducible representations of the rotation group present in 
this operator must destroy I^gs) separately. The scalar 
component is identically zero. The vector component is 



N 

E — 



z a - Zp 



[i(S a x Sp) + Sp]\^ GS ) =0 . (D3) 



Since I^gs) is a lso its own time- reverse it must be de- 
stroyed by the time-reverse of the vector operator, i.e. 



7 * _ 7 * 



[i(S a X,Sp)+ Sp] 



N 

E 

P^a 



Z 



Z a - Z, 3 



[i(S a xSp) + Sp] 



(D4) 



The difference of these is the trivial operator S a x S, and 
their sum is 2D a . 

We prove, now, the factorizability of Tins- In order to 
do so, we need the following identities: 



JV 



J2[i(S x S a ) + S] ■ D a 



In this Appendix we will prove the factorization for- 
mula, Eq. (p^) . In order to do so, we split the proof in 
two steps. First, we will show that the operator D a anni- 
hilate I^gs), then, we will prove the factorization equa- 
tion. Let us begin with the first proof. The operator 



N 



P^a 



^^Ks--(s a + h(si + h 

t z (3 £ * 



N 

J2iiS-(S a xD a ) + S-D a ] = -S-A , (D5) 



JV 

x + &\ ■ *S) + S] 



N - 1 



(S z a + x)] 



(Dl) 



^[N-1]S* , 



(D6) 
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N 



E 



Sf3 • S-y 



(z* - z*)(z a - z p ) 



N 



:S 2 



S a ■ Sp 

s ~^p\ z *- z e\ 2 



N 



■ (Sg X Sp) 



(D7) 



APPENDIX E: THE NORM OF THE STATES 

In this Appendix we provide the proof of the formula 
for the norm of the one-spinon and of the two-spinon 
eigenstates. 

Throughout all this section and the following one, we 
will make use of the scalar product between symmetric 
polynomials f(z\, . . . , Zm), defined as 



(f\9)= E f*(zi,...,z M )g(zi,...,z M ) 



(El) 



N 

' E 



(z a - z 1 ){z a - zp) 



, ^ z^z, \ { s a xS 1 )-Sp=A-S . (D8) 



By putting together the identity 

N N N 



EEE' 



i(S a x Sj) + gy]t ■ [i(S a x Sp) + Sp] 
(z* - z*)(z a - zp) 



N 



= Y^D a -D a + ls-K+l(N-l)S 2 (D9) 



and the identity 

N N N , 

EEE- 



i(S a x g 7 ) + gy] ■ [i(5 a x Sp) + gg] 
(z* -z*)(z a - zp) 



N N N 

-TTT - 



[Sj-Sp+iSj-(S a xSp)} 



3E 



S a • Sp N(N 2 + 5) 



\z a - zp\ 2 
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(D10) 



the proof is complete. 



Let us begin with the one-spinon eigenstates. The state 
for one spinon in coordinate space, ^t a (z\, . . . , zm), has 
been defined in Eq.@, in the odd-iV case (M = (N - 
l)/2) to be: 



,1/ 



M 



M 



.Z M 



)=Y[(z a -z J )Y[(z l -z J fl[, 



i<j 



(E2) 



is of the form $ Q x ^gs, where $ Q = Jl^C 2 " — z j). 
We will prove the formula for the norm of the state 
|$ m ), Eq.(|48|), by mathematical induction. In order to 
do so, let us define the symmetric operator: 



ei(zi, . . .,z M ) = Z\ + . .. + z M 



(E3) 



For any wavefunction of the form <!> x ^gs, where $ is a 
symmetric polynomial with degree less than N — 2M + 2, 
we have 



Hhs^gs = Egs^gs + ^(^") 2 *gs| \ 



E 



: 2 —. 

3 dz. 



^ Za - Zk 3zj 



and thus 



N -3^ d 
9 2^ 



3 dzi 



, (E4) 

j 



J 2-7T 

Unseizes ~ e x riHS^GS = ^(^) 2 



x * 



as 



E*f 



d N -3 



dz 



-ei 



(E5) 



Since ($|-D£ • D a \Q) is nonnegative for any wavefunc- 
tion |$), this provides an explicit demonstration that 
|\&gs) is the true ground state. The annihilation op- 
erators and their equivalence to TLhs when squared and 
summed were originally discovered by Shastry |Dfl . They 
are lattice versions of the Knizhnik-Zamolodchikov op- 
erators known from studies of the Calogero-Sutherland 
model, the 1-dimensional Bose gas with inverse-square 
repulsions pIMO. 



From the definition of $ Q one then obtains: 
M d d 



Eq.(E6) implies the following identity for \& r 

Eqs 



(E6) 



(J2tt 2 )/N' 



r][ei* m ] = [m(M - m) + M][ex%, 



1G 



+ [Af-(m-l)]* m _i 



(E7) 



(H is the scaled Hamiltonian: H = U H s I '(J2ir 2 /N 2 )). 
Then, the following identity chain is proved: 

(to - 1)(M - to + l)<* m _i|ei|* m ) 



:>,....ii;7/- v TO f^ /Ar2 ]|[ei^m]) 



(J2n 2 )/N' 



Y[M + l]Y[m +\]T[M -r, 
T[M + |]r[m + l]r[M - m + 1] 



2 Cm 



(E14) 



The constant Cm is expressed in terms of Wilson's inte- 
gral as 

M 

cm = (^/^.../^II( 1 --) 2 

2tti / zi / z M fj 1 : 2j 



= (M + to(M - TO))(* m _ 1 |e 1 |* m ) 

+ (M-m + l)(* m _i|* m _i) , 
thus, Eq.(|E^) implies the identity: 



= TV 



Af(2M)! 



(E15) 



(E8) Eqs.(E14 E15| ) complete the proof. 

Now we work out the formula for the two-spinon eigen- 
states. In this case M = N/2 — 1 and Eq.(E5) becomes: 



M -rn+1 



2(M-m+i) 



(E9) 



In order to determine a suitable induction relation, let us 
introduce the operator 



eAf(«ij • • • j Zm) = zi ■ . . . ■ z M 



Clearly 



e M (zi, ■ ■ .,ZM)eM{—> •■■> — ) = 1 
zi zm 



(E10) 



Since all the ^m's are products of the ground state factor, 
^GSi times a symmetric polynomial of degree less than 
2 in each variable, we have 



(# m -i|ei|tf m ) = ([e^tf m ]|ei|[e^tf m _i]) 



M- 



|ei|*Af-m+l) 



At this point, we use again Eq.(p9|) in order to write 
Eq.flElll) as 



(*M-m|ei|*A/-m+l) 



2(m-i) 



JT M-m\^ M-rn) 



,1/ 



(E16) 



where again we work with the Hamiltonian in scaled 
units. Eq.(p^) now becomes: 



M 







Mz„ - z z a — + Mzp - zl — 
az a p OZ0 



® a p , (E17) 



where the state has been defined in Eq. (|50|) . Hence, 
by letting e\ act onto the two-spinon eigenstate & mn , we 
obtain 



(EH) [H 



Eqs 3 
^ S /A . 2 ][ei$mn] - (A m „ + M - -)[ei$„ m ] 



(J2tt 2 )/N 2 







N 

ran \ ^ 



fc Q,/3=l 



(4) m W 



2[ei* a/3 ] 



2(m-i) 



Eq.(E12) closes the induction relation: 

QPmlttm) = (m-|)(M-TO+l) 

(* m _i|* TO _i) m(M-m+|) 

The formula generated by recursion is: 

(j-|)(M-j + l) 
j(M-m + i + |) 



(E12) 



(E13) 



(*m|* m )= n 



c 



j=l...m 



(Mz a -zl-^- + Mz,- Z j^-)^ 



(E18) 



which implies 



[H 



E, 



GS 



(J2n 2 )/N 2 -(^n + M+-)]e 1 ^ mn ] 



2aj^{[M-ro-fc+l]^&J ,+fc - 1 ' n - fc $ 



m-\-k— X-\-j,n— k— j 
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+ [M - n + k + 1] ^ b 3 



\ y Trn+k.n—k—l 



Qm+k+j.n-k-l-j} 



(E19) 



(See eqs.(56 8(]) for the definition of the coefficient 
af n ,bf n and eq.@ for the definition of X mn . 

L et us , now, take the scalar product of both sides of 
Eq.( E19 ) with & m -i. n . The result can be recast in the 
form 



($m-l,n|ei|$r. 



M-m+1 



($ m -l,n\$m-l,n) 2(M - m + h] 



(E20) 



On the other hand, by taking the scalar product of 
®m,n-i with both sides of Eq.(E19), we obtain: 



(A m , n _l — Xmn — M — — ) (& m ,n-l\ei\&mn) 



J2ar{[M- m -k + l} 



Hence, the proof is given by the following induction re- 
lations: 

($mn|*™n) 



($m,n-l|^m,7i-l) 



(n - |)(M -n+ |)(m - n + l) 2 
n(Af — n + 1)(to — n + | )(m — 71+5) 



($mJ$mn) 



($m-X,n|^m-l,n) 

(M — to + 1)(to — n + |)(m — n — \)m 
(M - to + |)(m - n) 2 (m + |) 

Eqs.( ^25p2^ ) imply 



(E25) 



(E26) 



mn I *mfi/ 



Ejm+fc-l.n-fe/j; i,*. 
Oj ^m.n-ll^m+fe-l+^n-fe-j 



(7 



, r[m - n + i]r[m - n + |] 



A/" 



T 2 [to - n + 1] 



+[M-n + fc+l]^ b ™+k,n-k-i 
j 

X ($m,n-l|$m+/!+j,ra-fc-j-l) j 

which implies the relation: 

($ m ,„-i|ei|$ mn ) 



($ m ,„_i|$ 



(E21) 



r[m + l]r[M - m + i] T[n + ±}T[M -n + 1] 
T[m + |]r[M - to + 1] I> + l]r[M - n + §] 



and the constant C M is now given by: 
, _ M (2M)\M+l 



(E27) 



(E28) 



m.n — 1/ 



(Af — n + |)(to - n + 1) 2 
2(to - n + |)(to - n + i)(M -?i + l) 



(E22) 



In order to complete the proof, we need two more iden- 
tities, which can be proved in the same way we did for 



Eq.(E13): 



($m-l,n|ei|$ m 
($mn|$mn) 



— \&M-n,M-m\ e l I $M-n,M-m+l) 



(m + i)(TO — 7i) 5 



2(m — n + i)m(m — n — i) 



(E23) 



($ m ,n-l|ei|$ T 



(^M-n,M-m|ei | $Af-ji+l,Af-m) 



($mn|$mn) 



2(n-|) 



(E24) 
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